Dynamics of the inner crust of neutron stars: hydrodynamics, elasticity and collective modes 
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We present calculations of the hydrodynamics of the inner crust of neutron stars, where a superfluid neutron 
liquid coexists with a lattice of neutron-rich nuclei. The long-wavelength collective oscillations are combi- 
nations of phonons in the lattice and phonons in the superfluid neutrons. Velocities of collective modes are 
calculated from information about effective nucleon-nucleon interactions derived from Lattimer and Swesty's 
microscopic calculations based on a compressible liquid drop picture of the atomic nuclei and the surrounding 
neutrons. 



I. INTRODUCTION 

Our focus in this paper is on the inner crust of a neutron 
star which, despite the fact that it occupies rather a small frac- 
tion of the total volume of the star, plays an important role 
in determining observable signals from neutron stars Q]. For 
example, heat from the interior of the star passes through this 
region on its way to the surface [2|, torsional oscillations of 
the crust have been proposed as a mechanism for quasiperi- 
odic oscillations seen in the afterglows of giant X-ray flares 
on neutron stars fl3}Q, and the properties of vortex lines in 
the neutron superfluid have been invoked as a source of glitch 
phenomena observed in pulsars. Collective modes of the in- 
ner crust have been proposed as a means of heat transport 
(6). A pioneering calculation of modes of the inner crust 
was performed by Epstein [7]. An important physical effect 
is "entrainment", the fact that not all neutrons participate in 
the motion of the superfluid, since some neutrons are effec- 
tively locked to the protons. This effect is familiar at lower 
densities in the outer crust, where there are no interstitial neu- 
trons and all neutrons are effectively locked to the protons. A 
general treatment of the hydrodynamics of a neutron super- 
fluid coupled to an elastic solid has been given in Ref. [8 |, a 
macroscopic treatment of entrainment has been given by [9|, 
a general formulation of the properties of crustal modes have 
been described in Ref. ifTOl . and a field-theoretical description 
of modes has been given in Ref. [ 1 1 1. Recently, Chamel, Page 
and Reddy have calculated frequencies of collective modes of 
the crust and have investigated some consequences of these 
modes [ 12]. 

In this paper, we consider hydrodynamics and collective 
modes of the inner crust. One of our main purposes here is to 
calculate velocities of long-wavelength collective modes. As 
shown in Ref. 1101 . key ingredients are derivatives of neutron 
and proton chemical potentials with respect to particle densi- 
ties. Calculating them is a challenge, since for the composite 
system of crystal lattice and superfluid neutrons it is necessary 
to take into account the equilibrium between neutrons within 
nuclei and those outside. In this paper we calculate the deriva- 
tives from results of Lattimer and Swesty lfT3ll for thermody- 
namic quantities; their microscopic calculations of dense mat- 
ter explicitly allow for the coexistence of nuclei and neutrons. 
Another important quantity is the neutron superfluid density, 



and we demonstrate its strong influence on mode frequencies, 
which has also been pointed out in Ref. lfT2l . 

To set the scene, we now give a few basic properties of 
crustal matter in neutron stars. It is a very good first ap- 
proximation to take the temperature to be zero, since it is low 
compared with the Fermi temperature of the electrons, jJ, e /k.B 
(where \i e is the electron Fermi energy and kg is the Boltz- 
mann constant) and with the melting temperature of the lattice 
for most of the life of a neutron star. The melting temperature 
T m of a Coulomb lattice of ions containing Z protons each is 
given by 
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where a t = [3/(4a:«;)] 1 / 3 = (9nZ/Ay^/k e (with k e = 
)'/ 3 being the electron Fermi wave number) is the ion 
sphere radius and tij, (j = n,p,e,i) are the densities of neu- 
trons, protons, electrons, and ions. The melting temperature 
will be affected by the presence of "dripped" neutrons in the 
inner crust, but we expect Eq. ([T]i to provide a good first es- 
timate even for this case. For an electrically neutral plasma 
consisting of a single species of ion and a rigid background of 
negative charge, the interaction parameter for melting, T m has 
been computed to be approximately 175 [ 1 1 and therefore 
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where a — e 2 / lie is the fine structure constant. This estimate 
shows that, given how rapidly neutron stars cool after forma- 
tion J2|, the ions are a crystalline solid for most of the life of 
the star. 

The remainder of this paper is organized as follows. In Sec. 
|II]a description of long-wavelength dynamics is given in terms 
of hydrodynamics and the theory of elasticity. Expressions for 
velocities of collective modes are derived in Sec. [ill] Deriva- 
tives of chemical potentials of neutrons and protons with re- 
spect to the densities of these components are evaluated in 



Sec. IV from Lattimer and Swesty's microscopic calculations 
of the properties of matter in the inner crust 1131 . There we 
also comment on the other important parameter determining 
mode velocities, the neutron superfluid density, and show that 
scattering from the ionic lattice could be reduced as a result of 
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ionic motion. We also give a simple expression for the shear 
elastic constant when the effect of screening by electrons is 
taken into account. In Sec. [V] we calculate velocities of col- 
lective modes and, in particular, show that they are strongly 
affected by the magnitude of the neutron superfluid density. 
Section[Vl]takes up a number of more detailed topics, and de- 
scribes outstanding problems and directions for future work. 

II. HYDRODYNAMIC FORMULATION 

The inner crust of a neutron star is a three-component sys- 
tem consisting of protons, neutrons and electrons. The basic 
long-wavelength modes are plasma oscillations of the elec- 
trons and two coupled modes in which the neutrons move in 
phase or out of phase with the protons and electrons, which are 
effectively locked together by the Coulomb interaction. Such 
modes may be described in terms of the neutron and proton 
densities, the ion velocity and a quantity specifying the motion 
of the superfluid neutrons. The proton, neutron and electron 
densities vary on length scales comparable to the lattice spac- 
ing and the quantities of interest for long-wavelength modes 
are coarse-grained averages of particle densities taken over 
distances large compared with the lattice spacing but small 
compared with the wavelength of the mode. The superfluid 
neutrons are characterized by an order parameter correspond- 
ing to l'(r) = (v/|(r)y/|(r)), where y CT (r) is the annihilation 
operator for a neutron with spin a and (. . .) denotes a thermal 
average. In the presence of the crystal lattice, ^(r) will gener- 
ally exhibit spatial structure on the scale of the lattice spacing. 
For a macroscopic description it is therefore convenient to in- 
troduce a phase 0(r), which is one half of a coarse-grained 
average of the phase of ^(r), the factor of one half being due 
to the fact that ^(r) corresponds to the wave function for a 
pair of neutrons with opposite spin. We then define the neu- 
tron superfluid velocity as 

vl = -ve. (3) 

m 

In this paper we shall extend the standard approach for deriv- 
ing two-fluid equations for liquid 4 He [14] at nonzero tem- 
perature and for mixtures of 3 He in liquid 4 He (TB |. We shall 
confine ourselves to long-wavelength phenomena and exclude 
short length-scale structures such as vortex lines. 

We begin by considering flow of neutrons when the ions 
are at rest. Phenomenologically, we can express the neutron 
current density to first order in the superfluid velocity as 

(j«)o = <«)o, W 

and the neutron momentum density is (g*)o = tnn s n \ s n . The 
quantity n s n is the neutron superfluid (number) density, which 
is related to the neutron superfluid (mass) density p;, as usually 
defined by the relationship p* = mn s n More generally, when 
nonlinear effects are taken into account, the neutron superfluid 
density n s n will be a function of v s n but in the linear theory, 
which is the subject of this paper, the dependence on may 
be neglected. 



Let us now perform a Galilean transformation to a frame in 
which the velocity of the protons (that of the ionic lattice) is 
x p . In this frame, the wave function of the neutrons acquires 
a phase factor exp(imv p • r/ti) at point r and therefore the su- 
perfluid velocity Q in the new frame is 

K) = K)o+v p . (5) 

The neutron current in the new frame is given by that in the 
old frame by the standard result for a Galilean transformation, 

j« = (jn)o+n„\ p =n s n (\ s n \) + n„\ p 

= «+»>„, (6) 

where n n n = n„ — n s n is the neutron normal (number) density. 
Equation |6]l has the standard form for the two-fluid model, the 
velocity of the lattice playing the role of the normal fluid ve- 
locity. Tunneling of protons between different nuclei is unim- 
portant because the proton chemical potential is well below 
the threshold for proton drip and therefore the current of pro- 
tons is given by 

jp = n p \ p . (7) 

In a two-fluid model, the velocity \ p corresponds to the veloc- 
ity of the normal component. 

We now give a hydrodynamic description of long- 
wavelength modes. We assume that matter is electrically neu- 
tral (n e = n p ). We shall further assume that the inertia of the 
electrons may be neglected, in which case the electrons re- 
spond instantaneously to changes in the proton density. This 
is a reasonable approximation, since the effective mass of an 
electron, p e /c, is small compared with the rest mass of a nu- 
cleon. The basic variables in a two-fluid description are then 
the densities of neutrons and of protons, the ion velocity, and 
9, the phase of the neutron condensate. We shall assume that, 
on the timescale of the modes, weak interaction processes that 
convert neutrons into protons may be neglected, and therefore 
the number of neutrons and the number of protons are each 
conserved. 

Conservation of the number of protons is expressed by the 
equation 

d t n p + V-(n p v p ) = 0, (8) 

and 

d,n n + V • j„ = d,n„ + V • «v„) + V • «v p ) = 0. (9) 

The total momentum density is g = g p + g„ and therefore the 
condition for momentum conservation is 

d,g + vn = o, (io) 

where II is the momentum flux density tensor. The final equa- 
tion is the Josephson relation, which gives the time evolution 
of the phase of the condensate. This reads 

hd t e+Hn = o, (ii) 

where \i n is the neutron chemical potential. 
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We now linearize these equations for small deviations of 
the system from a uniform state with no strain, a stationary 
lattice, zero neutron superfluid velocity and neutron chemical 
potential (jU„)o- The continuity equations <fsT> and (19} become 



d,n p - 

and that for neutrons by 



<? ( «„+n*V-v„+<V-v p = 0. 



(12) 



(13) 



In all cases, quantities not acted on by derivatives may be eval- 
uated in the uniform state. 

To identify contributions to the momentum flux density ten- 
sor it is convenient to consider the energy density of the sys- 
tem for small variations from the uniform state. The crystal 
lattice in the ground state is expected to be body centered cu- 
bic (bcc), but in a star one expects matter to be polycrystalline, 
rather than one single crystal. For wavelengths large com- 
pared with the characteristic size of crystallites, the medium 
behaves as if it were isotropic, in which case the deviation of 
the energy density E from its value Eq in the equilibrium state 
is given by 
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and L = —d 2 E / dn n duu. Changes in the proton density are 
related to the deformation vector by the relation 



n p (r + u) 



/(r + u,r)' 



where J is the Jacobian determinant, which for small strains is 
given by J ~ (1 + V • u). When linearized, Eq. ( 16 1 becomes 



Slip = -n p V ■ u = -ripUu. 



(17) 



It therefore follows that L = n p d 2 E / dn n dn p . The quantity 
K is the contribution to the bulk modulus due to lattice dis- 
tortions but without changes in the neutron density, and is 
given by K = n 2 d 2 E/dn 2 and S is the effective shear mod- 
ulus. The first two terms in Eq. ( 14 (are the usual expression 
for the elastic energy fl6l §4] of a solid, the third represents 
the interaction between density deviations of the protons (in 
nuclei) with those of the neutrons, while the final term is the 
self-interaction of neutron density deviations. First order con- 
tributions to the change in energy vanish because the initial 
state is taken to be in equilibrium. 

Changes in II come from two sources: one is an isotropic 
contribution 8Pdj k and the second is a shear contribution. The 



change in pressure is given by the Gibbs-Duhem relation, 
which at zero temperature is 



8P = n n dpL„ + n p 8pL p + n e 8pL e 
= n n Sn„+n p 8n epi 



(18) 
(19) 



where the latter expression holds for electrically neutral mat- 
ter and n ep = fi e + n p lf3T1 . Therefore the changes in II are 
given by 

8Yl ik = (n„Sli n +n p SHep)Sik-2S(uj k -^8 ik uii). (20) 

Equation ( [T0| for momentum conservation then becomes 
m(n p +rin)d t \ p + mn s n d,\ s n + n„V \l n + npWjX ep 

-2SVj{u ij - l -8 ij u ll ) = Q. (21) 



To linearize the Josephson relation (111 we introduce the 
variable 56 = 9 — (jU„)o t, which leads to the equation 



hd t S0 + Sn„ =0, 
or, taking the gradient of this equation, 
md,\ s n + V8fl„ = 0. 



(22) 



(23) 



On subtracting n s n times this equation from Eq. (21 1, one finds 



m(n p + n")d t v pi + n" n V ' i\i n + n p \ \\l ep 
-2SVj(ujj - ^8ijUu) = 0, 



m(n p + n n n )d t \ p + n"Vfi n + n p Vjj, ep 
4 

— SVV-u + SVx Vxu = 0. 
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(24) 



(16) F rom Eq- < fl4"l > one sees that 
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We show in the Appendix how the equations of motion (|8j, 
( 23 1 and ( |25] > may be derived from a variational principle. 



III. COLLECTIVE MODES 

The frequencies and eigenfunctions of long-wavelength 
modes may be determined by solving Eqs. ([12}, (13i, (23 1 
and p5j. We shall assume that physical quantities vary in 
space and time as e 1 ( q ' r_<i "), where q is the wave vector of the 
mode and co its frequency. Since we take the medium to be 
isotropic, there is no coupling between longitudinal and trans- 
verse modes. 

First, we consider transverse modes (q ■ u = 0). These in- 
volve no density changes and the frequencies may be obtained 
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directly from Eq. ( 25 i and one finds for the velocity v t = (o/q 



of the mode the expression 
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where p" = m(n p + n") is the total normal mass density. This 
result, which has previously been pointed out in Refs. iflOl 
[T21 . is the same as for transverse waves in a single-component 
solid, except that the mass density entering is less than the 
total mass density m(n p +n„) since only the normal neutron 
density participates in the motion and the neutron superfiuid 



remains stationary. 

We turn now to longitudinal modes (q x u = 0). In Eq. p5j ) 
the last term vanishes and it is convenient to eliminate V • u 



(27) there by the use of Eq. ( 17 1. This gives 



m(n p + n n n )d t \ p + n n n VpL n + n p VjX ep 
4 S 

+ V8n p = 0. 
3 «„ 



(28) 



Fourier transformation of Eqs. ( 1 2 1, ( 1 3 1, ( 23 1, and ( 28 1 yields 
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(29) 



On eliminating 8n n and 8n p from this equation, one finds 

pV-^ nn -4S/3 -S m 



p V - S s - 



0, (30) 



where 



= E pp n p n p 



+ E pn {nln[ 1 +n k n n<)+E, m {n n n ) 2 . (31) 



The upper indices k and / refer to the normal (n) and super- 
fluid {s) components and p s = mn s n is the neutron superfiuid 
mass density. For the protons, which are normal, n' p = n p and 

4 = 0. 

To understand the physics of the modes, it is convenient 
to find the frequencies of modes when the lattice is forced to 
remain stationary, and when the superfiuid velocity is forced 
to be zero. The velocity of the mode for a stationary lattice is 
given by 



vl = 



P 7 



(32) 



This mode corresponds to a phonon in the neutron superfiuid 
(the Bogoliubov-Anderson mode) and its velocity is given by 
the same expression as that for phonons in a Bose-Einstein 
condensate in an optical lattice [ 17 1. The velocity of the lattice 
modes when the superfiuid is stationary is given by 



S 



P" 



(33) 



The "force constant" for the mode depends on the energy due 
to changes in the densities proportional to the normal densities 
of the two species, together with an additional contribution 
due to the shearing of the lattice. The effective mass density 
is the total normal mass density. This mode corresponds to 



an oscillation of the lattice and entrained neutrons but no mo- 
tion of the superfiuid neutrons. In general, the two modes are 
coupled and their velocities v = co/q are given by 



(v 2 -v 2 n )(v 2 -vl)-vi =0, 



(34) 



where 
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The mode frequencies are therefore given by 
,,2 




or 



Since 
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(36) 



(37) 



(38) 



one sees that, if the shear modulus is neglected, one of the 
mode velocities will become imaginary if one of the condi- 
tions for thermodynamic stability to long-wavelength varia- 
tions of the neutron and proton densities, detEij > 0, is vio- 
lated. 

The way in which we have written the equation for the 
mode velocities differs from that which is used in Ref. Ifl2l . 
We have chosen to include all effects of entrainment, which 
gives rise to a current-current (vector) interaction, in the ve- 
locities v„ and v p . As a consequence only a scalar (density- 



density) coupling between these modes remains in Eq. (34i. 
In Ref. [12], corresponds to our v n but their lattice phonon 
mode is different from our mode with velocity v p and the cou- 
pling between their two modes has a vector character. 



5 



IV. EVALUATION OF PARAMETERS 



This condition is equivalent to 



The basic assumption made in the hydrodynamic treatment 
is that matter remains electrically neutral. This is a good ap- 
proximation provided the wave number is small compared 
with the Thomas-Fermi screening wave number of the elec- 
trons 



krF = \/4ne 2 dn e /dp e 



% 



(39) 



Here n e is the electron density, \i e the electron chemical po- 
tential, and k e is the electron Fermi momentum. The second 
expression in Eq. ( [39} holds for ultrarelativistic electrons. The 
condition for electrical neutrality at long wavelengths holds 
provided the frequency is much less than the electron plasma 
frequency 



(0„ 
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(40) 



For given proton and neutron densities, the mode frequen- 
cies depend on the parameters E nn ,E pp ,E np , the shear elastic 
constant, and the neutron superfiuid density n s n , which we now 
discuss in turn. 



A. E nn , E„ p , and E 



Determination of the En for the composite system of nu- 
clei immersed in a neutron gas requires information about the 
equilibrium state of this system. This may be determined from 
microscopic calculations. We have chosen to use the work of 
Lattimer and Swesty [13], which is based on Ref. If22"ll . One 
reason for doing this is that the model takes into account equi- 
librium between neutrons in nuclei and those outside and uses 
nucleon-nucleon interactions that have been fitted to masses 
of laboratory nuclei and the properties of pure neutron matter. 
A second reason is that second derivatives are evaluated ana- 
lytically without resort to finite differences, thereby reducing 
numerical noise. A third reason is that it is one of the mod- 
els of dense matter most widely used in simulations. A fourth 
is that extensive tables of physical quantities for the model 
are available. The particular variant of the model is relatively 
unimportant at subnuclear densities, and the one we used cor- 
responds to the parameters S v = 29.35 MeV, K = 370.6 MeV, 
n s = 0.155 fin" 3 , B = 16 MeV, m* = 0.911 , m* = 1.064m„ 
and S' v = ll .77 MeV and L = 3S' V . 

The general formalism described earlier applies for arbi- 
trary neutron excess but one expects that in the inner crust of 
neutron stars matter is close to being in equilibrium with re- 
spect to weak processes and consequently we describe results 
for this case. The condition for weak equilibrium is that the 
energy to add a neutron be the same as the energy to add a 
proton and an electron, i.e., 



H n = p p +p e = p ep . 



(41) 



P = p n -p p -p e = 



1 dE 

nd¥ e 



= 0. 



(42) 



where n = n p +n„ is the total baryon density and Y e = n p /n 
is the fractional proton concentration. While p — in the 
unperturbed state, it is nonzero during an oscillation because 
weak equilibrium cannot be established on the timescale of the 
modes. We took a temperature of 0.01 MeV, which is small 
compared with typical nuclear energy scales. 

On the Lattimer-Swesty website one can find values of 
dp/dn = d 2 E/dn 2 and dp/dY e . In beta equilibrium, 



dp_ 
dY e 



d 2 E 
l dndY e 



(43) 



To determine all the En we need a third quantity, and James 
Lattimer kindly prepared for us a table of dp/dY e , where 
p = ji„ — n p , from which one can calculate the quantity 
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(44) 



where the second expression holds in beta equilibrium. The 
calculations of Lattimer and Swesty are for electrically neutral 
matter and therefore from their results it is possible to deter- 
mine the sum of the proton and electron chemical potentials 
but not each chemical potential separately. In their approach, 
the energy is expressed as the sum of two contributions, the 
first being that of a gas of noninteracting electrons and the 
second being the rest, which includes the effects of nucleon- 
nucleon interactions and the Coulomb interaction. What is 
given as the electron chemical potential is the value for a free 
Fermi gas, 



Pe 



(Pic 2 - 



-mlc 4 



) l l 2 ^hc(S% 2 n e ) 



1/3 



(45) 



and all contributions from the Coulomb interaction are in- 
cluded in jip. More generally, if one wishes to allow for the 
electron and proton densities to be unequal, the electron chem- 
ical potential should include the contribution — e times the 
spatial average of the electrostatic potential. In the present pa- 
per p„ and p e enter only as their sum, and therefore it does not 
matter how Coulomb contributions are apportioned between 
the two chemical potentials. 

The calculations of Lattimer and Swesty give second 
derivatives of the energy density with respect to the variables 
n and Y e , but the quantities we need are derivatives with re- 
spect to «„ and n p . By changing variables, one finds that the 
two sets of quantities are related by the transformation 

dp/dn 
rC x dpjdY e 
dp/dY e 

(46) 

In Table I we exhibit results for physical properties of 
crustal matter. For the mass number A we have taken the 
number of nucleons within the volume of the nucleus plus the 
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number of surface neutrons, and w is the fractional volume 
occupied by nuclei. One noteworthy result is that E pp is close 
to the contribution djj, e /dn e from a uniform electron gas at es- 
sentially all densities. The Coulomb lattice contribution is of 
relative order Z 2 / 3 a, is negative, and is always small in mag- 
nitude compared with djj, e /dn e . In addition, contributions to 
Epp from nuclear forces are numerically small but they are 
positive and overwhelm the lattice contribution at higher den- 
sities. The effective neutron-proton interaction E np , which is 
attractive, is appreciable. 

B. Shear elastic constant 

In high density matter, there are a number of contributions 
to the elastic constants. The first is the Coulomb energy due 
to interactions between nuclei. This contribution was calcu- 
lated for a rigid background of electrons in Ref. lfl8l and it is 
of order «/Z 2 e 2 /a,. A second contribution is due to the non- 
uniformity of the electron density and this is negative and of 
order kf F af times the result for a rigid electron background. 
As we remarked earlier, we assume that the medium is effec- 
tively isotropic. To relate the effective shear elastic constant 
of the polycrystal to those of the single crystal is a difficult 
problem, and a number of suggestions have been made lfl9l . 
In an astrophysical context, Ogata and Ichimaru |20| did this 
by averaging elastic constants over all possible directions, a 
procedure which Baiko adopted in his recent work [21 1. For a 
static lattice of ions and a rigid background of electrons Ogata 
and Ichimaru find 

zV 

Seoul =0.1194/!,- (47) 

and the contribution due to inhomogeneity of the electron den- 
sity is well fit for Z < 40 by 

Screen = -0.010Z 2 / 3 [1 + (m e c/p e ) 2 } 1 / 2 5 C ° Ul (48) 

For applications to the inner crust of neutron stars, electrons 
may be treated as ultrarelativistic. The shear modulus has a 
negative contribution from lattice phonons and this has been 
calculated in Ref. [23 1 in the absence of dripped neutrons. We 
find that for T ~ 10 s K, the phonon contribution is about 2 
orders of magnitude smaller than the contributions for a static 
lattice and consequently we neglect it. The total shear elastic 
constant is given by 

S = Seoul + ^screen 

= 0.1194«, Z ^ (l-0.010Z 2/3 ) , (49) 
where in the last term we have taken the ultrarelativistic limit. 



C. Neutron superfluid density 

A key quantity is the neutron superfluid density, n s n . When 
nuclei occupy a small fraction of space, one would expect on 



T=0 r=10 S K T = 5xl0 S K T = 10 9 K 




FIG. 1: The factor exp(— W(2k° ut )) by which the scattering of the 
periodic lattice is reduced by lattice vibrations for four different tem- 
peratures. The quantity is calculated for the maximum momentum 
transfer 2k° ut of free neutrons at the Fermi surface for a number den- 
sity equal to n° ut and is therefore the minimum value of the factor. 



physical grounds that this would be close to the density of 
neutrons outside nuclei, which we denote by n out . However, 
in a recent series of papers, Chamel has carried out Hartree- 
Fock calculations of the band structure of neutrons in the inner 
crust of neutron stars which indicate that n s n could be as low as 
n out /10 or less at matter densities 1/10 to 1/3 of nuclear matter 
density [24|. There are a number of physical effects, such the 
effect of lattice vibration on the band structure, which we shall 
consider later in this subsection, and the influence of disorder 
and impurities that need to be investigated in order to reduce 
the uncertainly in n s n . To examine the sensitivity of collective 
mode properties to n s n we have carried out two sets of calcu- 
lations, one with n s n = « out and the other with «^ as a function 
of density given by the results of Chamel's calculations. 

The calculations in Ref. [24| were carried out for a rigid 
lattice of nuclei and one may ask whether the motion of ions 
could affect this. In the harmonic approximation for the lat- 
tice and for ions that have no internal degrees of freedom, the 
strength of the lattice potential corresponding to a lattice wave 
vector G is reduced by a factor exp[— W(G)], which is the 
square root of the Debye- Waller factor as conventionally de- 
fined l26l . For a cubic crystal, 

W(G)^GV), (50) 
o 

where (u 2 ) is the mean square displacement of an ion. At zero 
temperature, for a bec lattice of ions and a rigid background 
of electrons, one finds in the harmonic approximation that 




where ©k,£ is the frequency of a phonon with wave vector k 
and polarization e, and n^ £ = l/[exp(/ia)| 4 , £ /kgT) — 1] is the 
Bose distribution function. Thus one finds for zero tempera- 
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ture that (27) 



(« ;r=o : 



3 hu-i 
2 mjCOp 



(52) 



where m_i sa 2.79855 is the average of the inverse of the 
phonon frequency over the Brillouin zone, in units of 1 / co P i. 
Here the ion plasma frequency (Q p j is given by 



AnZ 2 e 2 ni 



(53) 



where m; is the mass of an ion. At the beginning of this sub- 
section we pointed to the uncertainty in the value of the neu- 
tron superfiuid density. This leads to a corresponding uncer- 
tainty in the effective mass of a nucleus, which is p"/«; = 
m p (n — n s n ) /rij. So as not to overestimate (u 2 ) we take it to be 
Am p . For nonzero temperature, one may estimate the thermal 
contribution to (u 2 ) by using the Debye model for the trans- 
verse modes, which are taken to have a constant velocity, and 
putting all longitudinal modes at a common frequency. One 
then finds 



h 



-U-\ 



6 T 

-— Dl (a t T pi /T) 



1 



1 



(54) 



a, expfaTpi/T) - 1 , 

where D\(x) — x^ 1 f^dyy/ (e y — 1). By fitting the model to 
the second moment of the phonon frequencies and to the av- 
erage of the logarithm of the phonon frequencies for a bcc 
crystal with a rigid background of negative charge, one finds 
a, = 0.426548 and a, = 0.88412 fl] p. 82]. The expression 
( f54| > is a good approximation for all temperatures below the 
melting temperature. In Fig. 1 we present numerical calcu- 
lation of e~ w ^ 2k " U ) for a number of temperatures. The mo- 
mentum transfer 2k° ut = 2(37f 2 /i° ut ) 1 ' 3 is the maximum value 
of the momentum transfer for neutrons on a spherical Fermi 
surface, and therefore this calculation gives a lower bound on 
the value of e~ w for other relevant momentum transfers. This 
calculation indicates that the Debye- Waller factor can be im- 
portant in reducing scattering from the periodic lattice and its 
effect on neutron band structure and the neutron superfiuid 
density should be investigated in greater detail. 

V. MODE VELOCITIES 



The transverse and longitudinal sound velocities (Eqs. (27 i 
and ( 36 1) are plotted in Fig 1 . Velocities of longitudinal modes 
are indicated by solid lines, and dashed lines show the veloc- 
ities of the uncoupled modes, v„ and v p , given by Eqs. (32i 
and (33 1 . As one sees, coupling between modes is relatively 



weak. With increasing normal neutron density, the velocities 
of both modes are reduced: a higher neutron normal density 
decreases the velocity of lattice phonon, while a decreasing 
neutron superfiuid density decreases the velocity of sound in 
the neutron superfiuid. 



o.i 
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FIG. 2: (Color online) Velocities of long-wavelength modes as a 
function of nucleon density. Velocities of uncoupled modes are given 
by v n and v p , the corresponding velocities of the coupled modes by 
v±, and the velocity of the transverse mode by v f . Curves are for 
n" n = n° ut while the crosses show results for n s n taken from Chamel's 
calculation 1241. 




FIG. 3: (Color online) Expanded version of the plot in Fig. 2 for 
densities at which the longitudinal mode velocities approach each 
other most closely for n s n = ;i out (see text). The dots show computed 
values, which have been joined by straight lines. 



In Fig. 2 we show a magnified plot of the longitudinal ve- 
locities for n s n = n° ut for a range of densities around 6x 10~ 4 
fm~ 3 where the velocities of the two modes are close to each 
other. The coupling between the two uncoupled modes, Eq. 

and E„ p terms can- 



( 35 1, passes through zero because the E„ 
ceL 



The calculations demonstrate that the reduction of the neu- 
tron superfiuid density predicted by Chamel would have a dra- 
matic effect on both longitudinal modes, whose velocities are 
reduced compared with the values for n s n = n out , as has also 
been pointed out by Chamel et al. [ 12 1. For the Chamel re- 
sults for n s n the velocities we find are in reasonable agreement 
with the calculations of Ref . |[T2l . 
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TABLE I: Microscopic parameters of crustal matter. The quantity x m is the proton fraction of matter in the interior of nuclei. 



n[fnr 3 ] 


« om [fm- 3 ] 


fl,,[MeV] 


y. = n p /n 




w 


A 


Z 


Em [MeVfm 3 ] 


E„ p [MeVfm 3 ] 


E pp [MeVfm 3 ] 


dp.Jdn e [MeVfm 3 ] 


2.512E-04 


1.138E-06 


2.729E+01 


3.557E-01 


3.573E-01 


1.380E-03 


106.0 


37.9 


3.937E+04 


-4.281E+04 


1.061E+05 


1.018E+05 


2.818E-04 


1.877E-05 


2.772E+01 


3.320E-01 


3.557E-01 


1.451E-03 


107.0 


38.1 


2.395E+04 


-2.236E+04 


7.081E+04 


9.870E+04 


3.162E-04 


4.239E-05 


2.806E+01 


3.070E-01 


3.544E-01 


1.510E-03 


107.8 


38.2 


1.906E+04 


-1.963E+04 


6.772E+04 


9.632E+04 


3.548E-04 


7.051E-05 


2.838E+01 


2.832E-01 


3.532E-01 


1.567E-03 


108.5 


38.3 


1.551E+04 


-1.800E+04 


6.685E+04 


9.414E+04 


3.981E-04 


1.032E-04 


2.870E+01 


2.610E-01 


3.521E-01 


1.626E-03 


109.3 


38.5 


1.274E+04 


-1.668E+04 


6.664E+04 


9.206E+04 


4.467E-04 


1.407E-04 


2.903E+01 


2.405E-01 


3.509E-01 


1.686E-03 


110.0 


38.6 


1.054E+04 


-1.548E+04 


6.667E+04 


9.002E+04 


5.012E-04 


1.836E-04 


2.936E+01 


2.218E-01 


3.496E-01 


1.750E-03 


110.8 


38.7 


8.757E+03 


-1.435E+04 


6.672E+04 


8.800E+04 


5.623E-04 


2.323E-04 


2.970E+01 


2.047E-01 


3.484E-01 


1.818E-03 


111.6 


38.9 


7.315E+03 


-1.327E+04 


6.672E+04 


8.598E+04 


6.310E-04 


2.877E-04 


3.006E+01 


1.891E-01 


3.470E-01 


1.891E-03 


112.5 


39.0 


6.141E+03 


-1.224E+04 


6.660E+04 


8.395E+04 


7.079E-04 


3.504E-04 


3.043E+01 


1.749E-01 


3.457E-01 


1.969E-03 


113.4 


39.2 


5.181E+03 


-1.127E+04 


6.634E+04 


8.190E+04 


7.943E-04 


4.213E-04 


3.083E+01 


1.620E-01 


3.442E-01 


2.054E-03 


114.4 


39.4 


4.394E+03 


-1.035E+04 


6.592E+04 


7.982E+04 


8.913E-04 


5.014E-04 


3.124E+01 


1.503E-01 


3.427E-01 


2.147E-03 


115.4 


39.5 


3.746E+03 


-9.491E+03 


6.534E+04 


7.772E+04 


1 .000E-03 


5.918E-04 


3.168E+01 


1.397E-01 


3.411E-01 


2.248E-03 


116.5 


39.7 


3.211E+03 


-8.688E+03 


6.461E+04 


7.558E+04 


1.122E-03 


6.938E-04 


3.214E+01 


1.300E-01 


3.394E-01 


2.359E-03 


117.7 


39.9 


2.766E+03 


-7.942E+03 


6.372E+04 


7.342E+04 


1.259E-03 


8.087E-04 


3.263E+01 


1.213E-01 


3.376E-01 


2.480E-03 


118.9 


40.1 


2.397E+03 


-7.253E+03 


6.268E+04 


7.123E+04 


1.413E-03 


9.382E-04 


3.315E+01 


1.133E-01 


3.357E-01 


2.614E-03 


120.2 


40.4 


2.087E+03 


-6.618E+03 


6.151E+04 


6.902E+04 


1.585E-03 


1.084E-03 


3.370E+01 


1.061E-01 


3.336E-01 


2.762E-03 


121.7 


40.6 


1.827E+03 


-6.036E+03 


6.021E+04 


6.678E+04 


1.778E-03 


1.248E-03 


3.428E+01 


9.957E-02 


3.315E-01 


2.926E-03 


123.2 


40.8 


1.607E+03 


-5.502E+03 


5.879E+04 


6.453E+04 


1.995E-03 


1.432E-03 


3.490E+01 


9.365E-02 


3.292E-01 


3.107E-03 


124.9 


41.1 


1.421E+03 


-5.015E+03 


5.727E+04 


6.225E+04 


2.239E-03 


1.639E-03 


3.556E+01 


8.827E-02 


3.268E-01 


3.309E-03 


126.7 


41.4 


1.262E+03 


-4.571E+03 


5.567E+04 


5.997E+04 


2.512E-03 


1.872E-03 


3.626E+01 


8.340E-02 


3.243E-01 


3.534E-03 


128.6 


41.7 


1.125E+03 


-4.166E+03 


5.398E+04 


5.768E+04 


2.818E-03 


2.134E-03 


3.700E+01 


7.899E-02 


3.216E-01 


3.785E-03 


130.7 


42.0 


1.007E+03 


-3.798E+03 


5.223E+04 


5.539E+04 


3.162E-03 


2.428E-03 


3.779E+01 


7.499E-02 


3.187E-01 


4.066E-03 


132.9 


42.4 


9.045E+02 


-3.464E+03 


5.043E+04 


5.311E+04 


3.548E-03 


2.758E-03 


3.862E+01 


7.136E-02 


3.156E-01 


4.382E-03 


135.3 


42.7 


8.150E+02 


-3.161E+03 


4.858E+04 


5.084E+04 


3.981E-03 


3.128E-03 


3.951E+01 


6.807E-02 


3.124E-01 


4.737E-03 


137.9 


43.1 


7.363E+02 


-2.886E+03 


4.671E+04 


4.859E+04 


4.467E-03 


3.544E-03 


4.045E+01 


6.509E-02 


3.090E-01 


5.137E-03 


140.8 


43.5 


6.667E+02 


-2.637E+03 


4.482E+04 


4.636E+04 


5.012E-03 


4.010E-03 


4.144E+01 


6.239E-02 


3.054E-01 


5.589E-03 


143.8 


43.9 


6.049E+02 


-2.412E+03 


4.292E+04 


4.417E+04 


5.623E-03 


4.533E-03 


4.249E+01 


5.994E-02 


3.016E-01 


6.101E-03 


147.1 


44.4 


5.496E+02 


-2.207E+03 


4.102E+04 


4.201E+04 


6.310E-03 


5.120E-03 


4.360E+01 


5.772E-02 


2.975E-01 


6.681E-03 


150.7 


44.8 


5.000E+02 


-2.021E+03 


3.914E+04 


3.990E+04 


7.079E-03 


5.777E-03 


4.477E+01 


5.571E-02 


2.933E-01 


7.340E-03 


154.6 


45.3 


4.553E+02 


-1.853E+03 


3.727E+04 


3.784E+04 


7.943E-03 


6.514E-03 


4.601E+01 


5.388E-02 


2.888E-01 


8.092E-03 


158.8 


45.9 


4.148E+02 


-1.700E+03 


3.544E+04 


3.583E+04 


8.913E-03 


7.340E-03 


4.731E+01 


5.221E-02 


2.841E-01 


8.949E-03 


163.4 


46.4 


3.780E+02 


-1.561E+03 


3.365E+04 


3.389E+04 


1 .000E-02 


8.266E-03 


4.868E+01 


5.069E-02 


2.791E-01 


9.931E-03 


168.4 


47.0 


3.445E+02 


-1.435E+03 


3.190E+04 


3.201E+04 


1.122E-02 


9.303E-03 


5.011E+01 


4.930E-02 


2.738E-01 


1.106E-02 


173.7 


47.6 


3.139E+02 


-1.320E+03 


3.020E+04 


3.020E+04 


1.259E-02 


1.047E-02 


5.162E+01 


4.802E-02 


2.683E-01 


1.235E-02 


179.5 


48.2 


2.860E+02 


-1.216E+03 


2.855E+04 


2.846E+04 


1.413E-02 


1.177E-02 


5.320E+01 


4.684E-02 


2.626E-01 


1.384E-02 


185.8 


48.8 


2.604E+02 


-1.120E+03 


2.697E+04 


2.680E+04 


1.585E-02 


1.323E-02 


5.485B+01 


4.575E-02 


2.565E-01 


1.556E-02 


192.7 


49.4 


2.370E+02 


-1.033E+03 


2.545E+04 


2.521E+04 


1.778E-02 


1.486E-02 


5.657E+01 


4.473E-02 


2.502E-01 


1.755E-02 


200.0 


50.0 


2.156E+02 


-9.524E+02 


2.399E+04 


2.370E+04 


1 .995E-02 


1.670E-02 


5.836E+01 


4.378E-02 


2.435E-01 


1.987E-02 


207.9 


50.6 


1.962E+02 


-8.789E+02 


2.260E+04 


2.227E+04 


2.239E-02 


1.875E-02 


6.023E+01 


4.289E-02 


2.366E-01 


2.257E-02 


216.4 


51.2 


1.785E+02 


-8.114E+02 


2.128E+04 


2.091E+04 


2.512E-02 


2.106E-02 


6.217E+01 


4.204E-02 


2.293E-01 


2.575E-02 


225.5 


51.7 


1.626E+02 


-7.493E+02 


2.003E+04 


1.962E+04 


2.818E-02 


2.364E-02 


6.419E+01 


4.124E-02 


2.216E-01 


2.950E-02 


235.0 


52.1 


1.484E+02 


-6.923E+02 


1.885E+04 


1.841E+04 


3.162E-02 


2.653E-02 


6.629E+01 


4.049E-02 


2.136E-01 


3.397E-02 


244.9 


52.3 


1.358E+02 


-6.399E+02 


1 .774E+04 


1.726E+04 


3.548E-02 


2.977E-02 


6.849E+01 


3.980E-02 


2.051E-01 


3.936E-02 


254.9 


52.3 


1.249E+02 


-5.919E+02 


1 .670E+04 


1.617E+04 


3.981E-02 


3.339E-02 


7.079E+01 


3.917E-02 


1.960E-01 


4.596E-02 


264.4 


51.8 


1.156E+02 


-5.481E+02 


1.573E+04 


1.513E+04 


4.467E-02 


3.743E-02 


7.322E+01 


3.863E-02 


1.862E-01 


5.421E-02 


272.6 


50.8 


1.080E+02 


-5.083E+02 


1 .483E+04 


1.414E+04 


5.012E-02 


4.193E-02 


7.581E+01 


3.820E-02 


1.756E-01 


6.480E-02 


277.4 


48.7 


1.021E+02 


-4.724E+02 


1 .399E+04 


1.320E+04 


5.623E-02 


4.693E-02 


7.857E+01 


3.791E-02 


1.639E-01 


7.884E-02 


275.3 


45.1 


9.788E+01 


-4.406E+02 


1.322E+04 


1 .228E+04 


6.310E-02 


5.244E-02 


8.155E+01 


3.777E-02 


1.507E-01 


9.828E-02 


259.8 


39.2 


9.532E+01 


-4.126E+02 


1.253E+04 


1.140E+04 


7.079E-02 


5.849E-02 


8.476E+01 


3.780E-02 


1.357E-01 


1.267E-01 


220.2 


29.9 


9.446E+01 


-3.883E+02 


1.190E+04 


1.056E+04 


7.943E-02 


6.510E-02 


8.819E+01 


3.795E-02 


1.183E-01 


1.712E-01 


145.6 


17.2 


9.562E+01 


-3.665E+02 


1.135E+04 


9.751E+03 


8.913E-02 


7.241E-02 


9.178E+01 


3.813E-02 


9.828E-02 


2.466E-01 


45.7 


4.5 


1.011E+02 


-3.429E+02 


1 .086E+04 


9.003E+03 



the ones traditionally used in condensed matter physics, and 
they have the advantage that the kinetic energy is a sum of a 
term proportional to u 2 and one proportional to (V0) 2 , with 
no cross terms (see Eq. ( |A4| >). This is similar to what is done 
in discussing stability of two-component systems, when it is 



VI. DISCUSSION AND CONCLUDING REMARKS 

In the general formalism, we used as the basic variables 
the lattice displacement, u, and the phase of the condensate, 
9, in addition to the neutron and proton densities. These are 



9 



convenient to work in terms of the density of one component 
and the chemical potential of the other one |29|. This ap- 
proach brings out clearly from the outset the fact that there 
is just one additional parameter needed to describe the ki- 
netic energy of the coupled superfiuid-lattice system com- 
pared with the uncoupled system. 

In calculating thermodynamic derivatives we used the cal- 
culations of Lattimer and Swesty, which treat Z as a contin- 
uous variable, which depends on density and proton fraction. 
Physically, it would be more realistic to keep Z fixed, since 
the time scales for changing Z, which requires either weak 
interactions or major rearrangements of the ionic lattice, are 
very long compared with periods of typical oscillations of the 
crust. However, in calculating derivatives of thermodynamic 
quantities we have allowed Z to assume its equilibrium value, 
but, since the variations of Z in the inner crust are rather small, 
we expect only small changes in mode velocities if Z is held 
fixed. 

In the calculations of transverse modes, we have included 
the effects of electron screening on the shear elastic constant, 
following the calculations of Baiko [21 1. This leads to a re- 
duction of the shear elastic constant by about 10% and of the 
velocity of transverse modes by about 5% for typical crustal 
conditions. 

Mode velocities are sensitive to the neutron superfluid den- 
sity, and in this connection we have pointed to the possible 
importance of phonons in reducing the strength of the peri- 
odic lattice potential experienced by neutrons. 

In this paper we have discussed the hydrodynamics and 
elasticity of the inner crust of neutron stars based on micro- 
scopic calculations of the properties of the matter lfT3ll that 
take into account in a thermodynamically consistent way the 
two-phase nature of the matter. One particularly important 
conclusion is that the effective neutron-proton interaction is 
appreciable. 

There are a number of directions for future work. One is to 
extend calculations to shorter wavelengths, where screening 
of the Coulomb interaction becomes less important. Another 
is to calculate damping of modes due to either Landau damp- 
ing, if electron mean free paths are long, or resistive losses, if 
the mean free path is short. It is also desirable to make calcu- 
lation of the effective nucleon-nucleon interactions based on 
state-of-the-art microscopic input, such as is obtained from 
chiral effective field theories of interactions between nucleons 
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Appendix A: Variational principle 

The equation of motion may be derived in a Lagrangian 
formulation with the Lagrangian density 



_S? = -/!„< 



I p « u 2„ 1<(V0) 2 ~„>-V0 
2 P 2m y ' 



-E(n n ,n p )-S[ u 



(Al) 



where the superfluid velocity is given by Eq. Q. For sim- 
plicity, we use units in which h = 1 in the Appendix. The 
independent variables are 9, n n , and the lattice displacement 
u, in terms of which changes in the proton density are given 



by Eq. (JT7 1. 

Alternatively, one may use a Hamiltonian formalism: 6 and 
— n n are conjugate variables, since the momentum conjugate 
to 9 is pg = dJ>?/dO = —n„ and the momentum density con- 
jugate to the lattice displacement is 



p n u-nye 



(A2) 



and therefore the Hamiltonian density is 

„2 , 



U 



2p n 



m 



K)_ 2 

p n 



(V0) 2 



P«u z 



~Pu • V6+E(n n ,n p ) +S Uj - ) (A3) 

(A4) 



It is instructive to compare the formalism for a neutron su- 
perfluid moving in a lattice with that for a homogeneous mix- 
ture of neutron and proton superfluids, as in the outer core of 
a neutron star. It is then natural to start from a generalized 
Ginzburg-Landau approach, with a Hamiltonian density 



GL 



1 
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(A5) 



where pjj is a matrix generalization of the usual superfluid 
stiffness coefficient or superfluid density for a one-component 
system (see, e.g., Ref. Il25l ). The current densities of the two 
components are then given by 



m 



(A6) 



To make contact with the results for the case of a lattice of 
protons, it is necessary to work in terms of the proton velocity, 
which may be defined in terms of the proton current density, 
which is given by 



Jp 



GL 



Ppp 



V0„ 



Ppn 



V0„. 



(A7) 
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Thus the variations of the phase of the neutrons produce a 
gauge term in the proton current. The mean proton velocity is 
therefore 



Pp P V0 p +p pn V0 n 



(A8) 



By comparison of Eqs. ( A10 1 and (|6| one sees that the expres- 
sions are equivalent if one makes the identifications 



(All) 



m 2 n 



Ppp Ppp 

Thus the neutron current density is given by 



/ Pnp \ V0„ p np 
J« = Pnn H -n p \ p . 

V Ppp I m Ppp 



(A9) 



(A10) 



and 



tlpPnp 
Ppp ' 



(A12) 
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